DMRG Study of the Striped Phase in the 2D t-J model 
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Using the density matrix renormalization group (DMRG), we study the 2D t-J model at a hole 
doping of a; = I on clusters as large as 19 x 8. We find a striped phase which is consistent with 
recent neutron scattering experiments. We find that bond-centered and site-centered stripes have 
nearly the same energy, suggesting that in the absence of pinning effects the domain walls can 
fiuctuate. 

PACS Numbers: 74.20.Mn, 71.10.Fd, 71.10.Pm 

In the low temperature tetragonal (LTT) phase of 
Lai.6-2;Ndo.4Sr2:Cu04, the tilt pattern of the CuOe octa- 
hedra form lines of displaced oxygens parallel to the Cu-0 
bond directions. These lines are rotated by 90° between 
adjacent layers. At a filling of a; = i, superconductivity 
is suppressed and neutron scattering studies reveal 
a striped domain wall ordering of holes and spins which 
is believed to be commensurately locked by the tilt dis- 
tortion of the lattice. One model for this striped order 
is illustrated in Fig. 1(a). Here the charge domain 
walls are shown running vertically and centered along 
the Cu-O-Cu legs, although the phase information re- 
quired to determine whether the domains should be leg 
centered or bond centered (centered between two legs) 
is not known. As shown, the domains are separated by 
four Cu-O-Cu spacings and for x = ^ contain one hole 
per two 4x1 domain wall unit cells. This latter feature 
is at odds with one-electron Hartree-Fock calculations ||] 
which predict a domain wall filling of one hole per domain 
wall unit cell. The spins in the regions between the walls 
are antiferromagnetically correlated with a tt phase shift 
across a domain wall. When x ^ ^, superconductivity is 
found to coexist with a weakened domain wall ordering, 
suggesting a close connection between the two. 

Here we present results of numerical density matrix 
renormalization group (DMRG) Q calculations for a t- 
J model with a hole doping x — ^. We find evidence 
for domain walls with tt phase-shifted antiferromagnetic 
regions separating the walls, and with a filling of one 
hole per two 4x1 domain wall unit cells. Kivclson and 
Emery Q have suggested that domain walls arise when 
phase separation of the holes into uniform hole-rich and 
hole-poor regions is frustrated by long-range Coulomb 
forces. The question of whether, in fact, the t-J model 
exhibits phase separation for the relevant physical val- 
ues of J/t and doping x remains controversial ||^,^ . Our 
present results show that long-range Coulomb forces are 
not necessary for the formation of domain walls. 

Depending on the dimensions and boundary conditions 
(BCs) of the cluster we study, the domain walls may be 
site-centered, as shown in Fig. 1(a), bond centered, or in 



FIG. 1. (a) Spin and hole structure suggested in Ref. [1] to 
account for neutron scattering experiments, (b) Hole density 
and spin moments for the central 8x8 region of a 16 x 8 
t-J system. The diameter of the gray holes is proportional 
to the hole density 1 — {ui}, and the length of the arrows is 
proportional to {Sf}, according to the scales shown. 

between. In contrast to Fig. 1(a), however, the site cen- 
tered domain walls have substantial hole densities over 
three rows of sites, rather than one. Previous attempts to 
understand the charge degrees of freedom of the striped 
phase have focused on one dimensional models . We 
consider another approach, in which coupled ladders are 
used to model the 2D system. In particular, in order to 
understand bond-centered stripes, we consider an array 
of two-leg ladders which are coupled antiferromagneti- 
cally via a mean field. We find that the 7r-phase-shifted 
magnetic order of the bond-centered striped phase can be 
understood within this mean field picture. Similarly, the 
magnetic order of the site-centered striped phase can be 
understood in terms of doped three-leg ladders coupled 
antiferromagnetically to undoped chains. 

The t-J Hamiltonian in the subspace of no doubly oc- 
cupied sites is given by 
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operator which creates (destroys) an electron at site i 
with spin s. The near-neighbor hopping interaction is t 
and the near-neighbor exchange interaction is J. We refer 
to the Cu-Cu lattice spacing as a and measure energies 
in units of t. We consider only J/t — 0.35 here. 

We present results here for L x 8 clusters, with L as 
large as 19. As first discussed by Liang and Pang the 
truncation errors in a DMRG calculation typically rise 
exponentially with the width of a large two-dimensional 
system (while only linearly with the length). However, 
the errors also tend to fall exponentially with the number 
of states kept per block. Consequently, while studies of 
doped L X 8 clusters are quite difficult, by keeping from 
1000-2000 states per block, we can obtain useful results, 
with truncation errors of 0.0002 - 0.0001. We are able to 
keep this many states because of recent improvements in 
the DMRG finite-system algorithm . 

The nature of the ground state of the 2D t-J systems 
causes additional numerical difficulties. Rather than an 
approximately homogeneous phase, we find that the sys- 
tem tends to have inhomogeneous charge and spin dis- 
tributions (such as domain walls), which can be pinned 
by the open BCs usually used in DMRG. Usually more 
than one such low-energy configuration is possible: for 
example, one could have horizontal as opposed to verti- 
cal stripes. A DMRG calculation involves sweeps through 
the sites of the lattice, and the energy of the approximate 
DMRG ground state of the system is decreased mostly 
through "local" improvements of the wavefunction. We 
find that in a large 2D system DMRG is usually unable 
to tunnel between two substantially different low-energy 
configurations. Even when a low energy tunneling path 
exists between two very different configurations, the cal- 
culation may move along the path slowly. To deal with 
these difhculties, we usually perform several simulations 
for each system. These systems differ in the charge and 
spin configurations in the first few sweeps. Later sweeps 
drive the system to a local energy minimum. One can 
then compare the total energy of different simulations to 
find which configuration is the ground state. The charge 
and spin configurations can be controlled in two ways: 1) 
by adjusting the total quantum numbers of the system 
at each step as the lattice is first built up from a few 
sites; and 2) by applying local chemical potentials and 
magnetic fields for the first few sweeps. Unfortunately, it 
is possible to miss the true ground state configuration if 
it is substantially different from what one expects. How- 
ever, unlike an ordinary variational calculation, only the 
crudest overall features of the wavefunction, such as the 
general location of the domain walls, are specified in the 
initial sweeps. These various runs can give substantial 
insight into what sorts of low-energy configurations can 
possibly occur under slightly different BCs or small per- 
turbations to the Hamiltonian. 

Figure 1(b) shows the charge and spin density in the 
ground state for the central 8x8 section of a 16 x 8 



system with J/t = 0.35 and 16 holes, corresponding to 
a filling x = Periodic BCs were used in the y direc- 
tion, and open BCs in the x direction. Along the left 
and right edges of the system a small staggered magnetic 
field of O.lt was applied. The BCs and the edge stag- 
gered field serve to orient and pin the domain walls in 
the configuration shown. In an LTT phase, the domain 
walls are oriented, and possibly pinned, by the lattice dis- 
tortion. The staggered edge field further acts to pick a 
direction for the spin order, which allows direct measure- 
ment of the spin configurations and reduces truncation 
errors in the DMRG calculation. Previous to this calcu- 
lation, dozens of simulations were performed, mostly on 
8x8 clusters, to find the nature of the ground state and 
the effect of various BCs. Included were several initial 
conditions corresponding to phase separation, with the 
hole cluster either on the edge or in the center of the 
system. These phase-separated configurations were un- 
stable, with the hole cluster tending to split or lengthen 
into domain walls. A single eight-hole vertical domain 
wall was also unstable, even when initial conditions and 
boundary staggered magnetic fields favored one. Objects 
resembling diagonal domain walls have been observed in 
three chain and four chain ladders ||l^,|l^, but attempts 
to stabilize a diagonal domain wall on an 8 x 8 system 
instead yielded a bent domain wall with the central part 
aligned in the (1,0) direction. Periodic BCs in the y di- 
rection tend to favor vertical domain walls; open BCs in 
the direction of the stripes tend to suppress them. 

In the simulation shown in Fig 1(b), eleven sweeps 
were performed, and in the final sweep 1400 states were 
kept. A local chemical potential was applied to confine 
the holes to the width-two stripes shown for the first six 
sweeps, and then removed. No initial magnetic field was 
needed away from the left and right edges to orient the 
TT-shifted antiferromagnetic domains as shown. 

Figure 2 shows the domain wall structure in a different 
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FIG. 2. Average hole density nh(ix) (solid circles) and 
spin structure function S-wiix) (open squares) for the 16 x 8 
system of Fig. 1(b). 
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way. With the sohd circles we show the local hole density 
nh{i) = 1— < cl^ci^i + cl^cei > as a function of the x- 
coordinate l^- The bond-centered nature of these domain 
walls is clearly evident. To show the spin structure, we 
define 

SM = {-ly^+h (2) 

With the open squares, we show STr{lx)- The period-8 
spin structure is clearly evident. 

The boundary conditions have a strong effect on the 
structure of the domain walls which appear. Bond- 
centered domain walls tend to form one lattice spacing 
away from an open boundary. This initially led us to be- 
lieve that sitc-ccntcred domain walls were not stable, but 
subsequent simulations showed that site centered walls 
occur also. In Fig. 3 we show the local hole density and 
spin structure function S-,r{lx) for a 19 x 8 cluster with 20 
holes. The same BCs and edge magnetic field as for the 
system shown in Fig. 2 were applied. A system such as 
this with an odd number of domain walls and open BCs 
in the x direction is forced by symmetry under reflection 
about a vertical line to have a site-centered domain wall 
in the center if is odd, and a bond-centered wall if 
Lx is even. In the calculation shown, reflection symme- 
try is used explicitly, which ensures that a site centered 
domain wall appears in the center. Note that the second 
and fourth domain walls, which are not so constrained by 
geometrical effects, are more site-centered than bond cen- 
tered. We have compared the local energies averaged over 
4x4 regions covering site-centered and bond-centered 
walls; the difference in energy per site between these was 
within our numerical errors for local energies, with both 
giving E/N w 0.62t ± O.Olt. In addition to bond cen- 
tered and site centered walls, asymmetrical walls can oc- 
cur. The close energy differences between these different 
types of walls suggests that a large 2D t-J system at 
a; = 1/8 might have fluctuating domain walls. 
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FIG. 3. Average hole density nh(£x) (solid circles) and 
spin structure function Stv{£x) (open squares) for a 19 x 8 
system. 



The TT-phase-shiftcd antifcrromagnctic regions reduce 
the energy for transverse hopping of holes within a do- 
main wall. To understand in more detail the bond- 
centered striped structure, we consider a model of antifcr- 
romagnetically coupled two-leg ladders. Ladders doped 
with X = 0.25 are alternated with undoped ladders, and 
no hopping is allowed between ladders. Ladders are 
exchange-coupled via a mean field, which is staggered 
along a ladder, but which may or may not have a tt 
phase shift across a doped ladder. The properties of a 
single ladder are calculated with DMRG, with a static 
magnetic field with wavevector (7r,7r) or (0,7r). In Fig. 4 
we show the magnetic response KiS'^)! to an applied field 
with magnitude h. As expected, an undoped ladder has 
a much greater response at the Neel wavevector (7r,7r). A 
doped ladder, in contrast, shows a substantially greater 
response at (0,7r). Hence the mean field treatment shows 
the TT phase shift seen in the 2D calculations. The mean 
field self-consistency conditions are 

hd,u = J\{S,)u,d\ (3) 

where u and d stand for doped and undoped ladders. 
From the results shown in Fig. 3, we find |(S'2)t,| ~ 0.32 
and \{Sz)d\ = 0.15. The results from the 16 x 8 system, 
in contrast, are \{Sz)u\ = 0.29 and |(S'z)d| = 0.13. As 
one expects, the mean field treatment overestimates the 
magnetic order. (We expect that correction for trunca- 
tion errors and finite size effects would further decrease 
the DMRG results.) The energy of this mean-field striped 
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FIG. 4. Magnetization per site \{Sz)\ induced by an ap- 
plied magnetic field h at wavevoctors (tt, tt) and (0, tt) on a 
2 x 32 ladder, (a) An undoped ladder, (b) A ladder with 
doping X = 0.25. 
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phase (with a; = and x = 0.25), including the exchange 
couphng between ladders, is about 2% higher than the 
energy of an array of uncoupled ladders at uniform den- 
sity {x = 0.125). Hence the mean field approach does not 
predict the charge ordering of the striped phase. Nev- 
ertheless, these results suggest that coupled ladders are 
natural starting points for understanding striped phases. 

A similar mean field treatment can be made for site- 
centered domain walls, coupling doped three-leg ladders 
with undoped single chains. This also yields 7r-phase- 
shifted antiferromagnetism, with reasonable magnitudes 
for KS'z)]. We will present these results elsewhere. 

According to the Maxwell construction, the occurence 
of phase separation results in a linear dependence of 
the energy on the filling over a range of fillings. While 
there has been disagreement in previous studies about 
whether phase separation occurs in the low-doping region 
for J /t ~ 0.3 — 0.5, it is clear that the curvature in the 
energy versus filling curve is small [Q,^. The possibility 
of a striped phase, which has generally not been consid- 
ered in these studies, makes the analysis more difficult. 
In particular, if we assume that at low doping, holes go 
into a single domain wall, then it would appear from the 
energy alone that one has phase separation up to a filling 
of ~ N~^^^, where N is the number of sites in the system. 
At higher dopings, an array of weakly repulsive, widely 
spaced domain walls would show a nearly linear energy 
versus filling dependence. In fact, we have observed ex- 
actly this scenario on a 12 x 6 system at J/t = 0.5, with 
periodic BCs in the y direction and open BCs in the x 
direction. Two holes bind into a pair, with a binding 
energy of 0.26i ± O.Oli. Two pairs bind into one vertical 
domain wall, with binding energy 0.10i±0.03i. These do- 
main walls are very similar to the ones seen in the L x 8 
systems, with either site-centered or bond-centered walls 
possible, but with a linear hole density greater by a factor 
of 4/3. Eight holes form two widely spaced domain walls; 
twelve holes form three domain walls. We will report on 
these results in more detail elsewhere. 

Why is superconductivity suppressed specifically at 
X = 1/8 in Lai.6-a;Ndo.4Sr^Cu04? More work is needed 
to answer this question, but we can make some general 
statements. First, viewing the stripes as coupled two- 
leg ladders, there does not appear to be any anomalous 
feature in an isolated two-leg ladder at a; = 0.25, such 
as a charge gap, which would suggest that domain walls 
must occur with exactly this doping. In any case, the 
charge on the domain walls spills out onto the rows of 
adjacent sites; our bond-centered domain walls have a 
hole density of x w 0.18 on the walls and x « 0.07 on 
the adjacent sites. This suggests, along with our results 
on the 12 X 6 system, that domain walls can occur with 
a range of dopings. 

The period of charge density wave (CDW) correlations, 
however, is insensitive to the broadening of the walls. We 
have observed this effect on width-four models of domain- 



walls. The CDW can be viewed as a one dimensional line 
of hole pairs. The pairs extend beyond the two-site width 
of the wall, but the period is set by the one-dimensional 
hole-density. This period is 4a at x = 0.25, the same as 
the transverse period of the stripes. In the LTT phase, 
the CuOg tilt structure causes the domain walls to be 
perpendicular in adjacent planes [Q. This suggests that a 
coupling between planes, such as through an electrostatic 
potential [|j or through a lattice distortion, could induce 
a static CDW order along the domain walls. This CDW 
order would tend to suppress superconductivity. 
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